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^ : Abstract 

> 

oo 
O 
(N 

Qv^ ' The partition function can be expressed as a power series in (3 (inverse temper- 

ature) . The leading term in the high temperature expression of the partition 

o ' function and effective potential is calculated exactly. Physical quantities like 

y^ ' the mean square separation can also be exactly determined in the high tem- 



perature limit. 

1. INTRODUCTION 

The study of string theory at finite temperature has received renewed attention recently 
P^r?| . In a recent paper one of us attempted to elucidate the nature of the Hagedorn 



transition |[1[] using the matrix model and found similarities with the deconfinement transition 
in gauge theories. This was also investigated in a subsequent paper using the AdS/CFT 
correspondence 0] . It is clear that recent developments in non-perturbative string theory or 
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M-theory |]T8|-Q have some bearing on our understanding of the high temperature behavior 
of strings. Furthermore, study of high temperature behavior of a system is often a useful 
probe of the system because many calculations simplify at high temperature. We can thus 
hope to learn something about M-theory from its high temperature behavior. With this 
motivation, in this paper, we attempt to study the high temperature behavior of a simple 
but non-trivial system - the system of two D-0-branes. This is essentially the BFSS matrix 



model |T^ with N=2. This is not big enough to describe M-theory. In particular it would 
not include for instances processes involving pair-production of D-0 brane - anti-D-0 brane. 
Nevertheless it is already complicated enough. In particular the nature of the threshold and 



other bound states that have been studied [p^ -p7[| are not fully understood. Furthermore 



we should keep in mind that while the matrix model reproduces string theory at short 
distances, the fact that it also does so at long distances seems to be entirely due to the 
super-symmetric non-renormalization theorems. At finite temperature supersymmetry is 
broken and perhaps we should not expect this. For all these reasons the study of matrix 
models at high temperature is worthwhile. 

A related model of D-instantons, the IKKT matrix model POl, which is 0-1-0 dimensional 



has been solved exactly for N=2 |^. The D-0-brane action that we are interested in, is 
a quantum mechanical one (i.e. 0-1-1 dimensional). However after compactifying the Eu- 
clideanised time, if one takes the high temperature limit, it reduces to a 0-1-0 dimensional 
model. There is thus a hope of solving this model order by order in (3 but to all orders in 
Qs using the same techniques as [§. One can then calculate physical quantities such as the 
mean square separation of the D-0-branes - a measure of the size of the bound states. This 
is what is attempted in this paper. We obtain the leading behaviour in j3. We can also 
estimate, the corrections to the leading result. The noteworthy feature being that each term 
is exact in its dependence on the string coupling constant . 

This paper is organized as follows. In section ^ we set up the problem and in section |^ 
we describe briefly the solution. The last section contains a summary of the result and some 
conclusions. 



2. THE ACTION 



The (0 + 1) dimensional BFSS lagrangian is 

9 _ 7 _ . . .1 

(2.0.1) 



L = — -tr 
2gls 



xix^ + 2tee - ^[x^x'^]2 - p^,[e,x^] - ^[x°,x^]x' 



where i = l,...,9;/i = 0, ...,9 ; X^ and 6 are N x N hermitian matrices. X^ is 10 a 
dimensional vector and ^ is 16 component Major ana- Wyel spinor in 10 dimensional super- 
Yang Mills theory. For X = 2, we can re- write X'^ and 6 as 

X^ = J2 l^^XJi, e = J: ^cr-Oa (2.0.2) 

a=l ^ a=\ ^ 

where, a" are the hermitian Pauli matrices and X^, Qa are real fields in terms of which , the 
Lagrangian ( p.O.lQ is 

1 f 3 . . 3 . 2? 3 

^ahc -| 3 1 3 I 

+ 7«"-^a-^6-^c +^ E Xa^b^a^b ~ l^M E ^a^b^b^l \ (2.0.3) 

If we Euclideanize and compactify time on a circle of circumference (3 the action becomes 

S = i I Ldt (2.0.4) 

Jo 

and the fields satisfy the boundary conditions 

x^(o) = x^ip), e{o) = -e{f3) 

Considering these boundary condition, we can expand the fields X^, 6 in modes as 

oo oo 

w — ^ 2Trin t ^ — ^ ZTrir > 

mt) = E Kne~\ Oa{t) = Y: Oa,re — ' (2.0.5) 



n=— oo r=— oo 



where, n is an integer and r is a half-integer. 

So, in terms of these modes, the action reduces to 



''s r,s,i=-oo, ^ *sP ri+Z+p=0 

1 OO 1 oo 

"*" 2/4 -^ -^a,n-^b,m-^a,l-^b,p ^ ^ Z^ -^a,n-^b,m-^b,l-^a,p ( (2.0.6) 

n+'m+l+p=0 n+»n+/+p=0 

n, m, I, p are integers and r, s are half-integers. 

Here, X^^, /3 have the dimension of length (L) and 9 has the dimension La . We scale X^^, /? 

-1 
with a factor of 7-, and 9 with a factor /s^ to make them dimensionless, which is equivalent 

to replacing all the Ig in the action by 1. 

The first and the second terms in the action give the masses of the modes of the vector and 

the spinor fields respectively, namely, ^ and %^. 

The partition function with this action is 

"e-'^ (2.0.7) 

In the next section we will try to calculate this . 

3. THE PARTITION FUNCTION AND THE EFFECTIVE POTENTIAL. 

3.1. Pfaffian 

The fermionic terms in the action are of the form {9a^rlolp,9b,sXcn)n+r+s=o and 9a^r9b,-r- In 
(3 —>■ limit, the first term in the action (p.0.6|) is dominant and (X^„)„^o is of the order 
^/J3. So, in large temperature limit the first term contributes to the partition function in 
the leading order of P only for n = {9a,rlolfi9b-rXc^o) ■ With these terms the action takes 
the form 

i(3 °° ( Anr Anr . .1 

•S*/ = -— Y] < —^9a,r9a-r ^9a-r9a,r — 2z6'a,r7o7At^6.-rX^O^'' '^ ~ 226'a ^r7o7^6'fc^X^of " '^ f 

Now, we will try to find out the pfaffian 



3 
a=l 



iSf 



(3.1.1) 



for this action 



CO „ 3 
r=Q-' a=l 



-^ i0a,r (^5^' - ^X,^oe"'^7o7.) 0,,. 



27rr 



<5"'' + ^X,V«^^7o7J^;,,r 



(3.1.2) 



We rotate X^g by a Lorentz transformation so that only X^q, X^q and X^q are nonzero. 
We take the representation of the Gamma matrices, in which 



7o = io-2 ® Is, 



7i = o^s ® Is, 



72 = -o"! ® Is 



(3.1.3) 



With this choice of representation we can write 



5/ = ^ E (^a,.f^5"^-^e'^^'=(x°l2 + X>i+X>3))®M, 



P 



b,—r 



(3.1.4) 



So, the pfafhan will be 

3 

n dO'^,rd0b,-r^^P 



'f 



afi=l 



^-{e'i^-^5^'-^e'^'^ 



(^°oi2 + xl^ai + xf 0O-3J J e[_.^. 



.9 I - V /3 
where 9'^ ^ and ^^ _^ are the spinors in three dimension, and has two components. 
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y r>0 



]J 16r H r— X^gX^o H — ((X^gX^o) - 4(X^gX^gX^gX^g) 



/3^ 



"4(X°gX°gX^gX^^g)j + — (^eafece'''''^X^ gXj'^gX, 



(3.1.5) 



We can see that the above expression has 5*0(3) symmetry in spinor indices, and 5*0(2, 1) 
symmetry in the vector indices. The 16r^ term gives the free fermionic contribution. Note 
that it is temperature independent as the Hamiltonian is identically zero for free fermions 
in + 1 dimensions. 

3.2. Free Bosonic sector 



After doing the fermionic integral, the partition function is 



Z=IlJd''X^^^ 



a=l 



exp 



916 24 / o3 (32 

TT 4r^ + — e ;, f^"''^X^' Xlf X^ + ——X^ X^ 
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E 



nj^O 
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rnie' 
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n+l+p=0 



I -■- V^ yn yv y^l yv ^ \^ yfl yv y^i yu 

n,m,l,p^ — QO, n,m,Z,p— — oo, 



(3.2.1) 



In infinite temperature limit ie (3 ^ 0, the first term dominates. To see the comparitive 
j3 dependence of the other terms; 

1) we set {Xji^jn^o ^ VPX^,n 

2) keep the terms contributing to the leading order of the partition function in /? ^ limit 
(this is justified in Appendix ^. 

3) transform back v^X^„ -^ (X^„)„^o 

and we can write the partition function up to a numerical factor 

1 



■'boson 



D \ri iD \ri 



— rf^x;,rf^x:„„exp 
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n=— oo;n^0 



(3' 



d^Xo Y[ d^Os exp 



— { -{Xafi.Xafi){Xbfl.Xb^Q) — -{Xafi.Xbfl){Xb^Q.Xafi] 



(3.2.2) 



Thus Zfjoson = /e ^"^ = ZfreeZo wherc the first part of the partition function Zjree is just 
the free bosonic particle partition function (per unit volume), which is 



3 9 



n a=l 1=1 



Zfree = I] IT IT / d"" XU"" X^,^ exp 



Y^ ^^' "' (VI Vt 



n^ — CX3; 



P' 






54 



Now, using E^=i n-' = ({s) and fX{s) = - E^=i ^"^ log n and C(0) = 1, C'(0) = -^ log(27r) 



ds 

(where ({s) is the Riemann Zeta function), we get 

27 
/ I \ 

^ free 



1 ' 

J~9. 



(3.2.3) 



Note that the free fermionic contribution contribution was discussed in the previous section. 
Zq is calculated below. 



3.3. Leading Interaction Term. 

Now we try to calculate the effect of the interactions. 
As argued earlier the leading (3 dependence is given by the zero modes, so in the first 
approximation we drop the terms in the action involving the higher modes. Thus we get, 

c r^ / \A/i yu -yP- vi^ vl^ yi^ yfJ' yv \ 

^P f/AA \2iv \2 , rv \2tv \2 I /-v \2/v \2 



~ (-^i,o--^2,o) — (-^2,0-^3,0) — (-^3,o--^i,o) I (3.3.1) 

We would like to first calculate the leading order contribution to Z that is 

Zo = |rfXV-'^" (3.3.2) 

Now, consider the parametrisation 

^1,0 = (a^i,^i), ^2,0 = (3^2,^2), X3,o = (^0) (3.3.3) 

If we had considered the original action with all the modes the action would not be Lorentz 
invariant, for example the terms X°;X^^X*„ and X*„X* _„ in the original action are not 
Lorentz invariant. However, expression (|3.3.1D has Lorentz invariance, and we are justified 



in using this parametrisation in order to evaluate ( p.3.2|) 



Under this parametrisation, the action (|3.3.1|) takes the form 

Sq = — |r?r2 sin^a + xfrl + r^Xn + r^P + l'^r\ -\ — 2xi.X2rir2 cosaj- (3.3.4) 

The partition function is 

Zo= ! d''Xl,d''Xl,d''Xl,e-^'- (3.3.5) 



At this stage, we can find the temperature dependence of the partition function and the 
mean square separation of two D-0 branes from simple scaling argument. As we have seen 
the leading order i.e. the zero mode contribution of the partition function comes from the 



term -^^[X^^ , X'^]'^ in the Lagrangian (|3.3.1|) . And the above parametrisation we have used 
here is also vahd for SU{N) matrix model, only the functional form of ^[X^, X'^]'^ in terms 
of this parametrisation will be different for different A^, but in each case the function will 
be homogeneous in /, Tj, Xj where < z < A^^ — 2 and of order 4. So, in each case we need to 
scale these variables hj P~^gi to scale out the /? from the exponent. And the temperature 
dependence of (/^) will be I3~^g2 in the leading order. Under the scaling above the measure 

15 15 

in Zq will pick up a (3~~ g~ factor for SU{2) , which comes from (3 x 10) X^^. In general 
for SU{N) in D dimension there will be D{N'^ — 1) X^„ in the measure. So, the partition 

g(JV^-l) £)(iV^-l) 

function Zq has temperature dependence (3 4, g 4, , And Zfree will be proportional 

to (/5(7)(^~i)(^'-i). 

Now we evaluate the partition function for this action. 



Zo = f dldQ^^'^l^ fdxidx2 f dridr2dQf^dn'^^darl4sm'^ ae-'^ 



(3.3.6) 



Using J dQ^^^ = 2*^" ^^tt, and after integration over xi , X2 and r2 the partition function is 

2^'^g^TT'^T{A) [^ ,„g f^ , 7 [ 13 



Zo 



/?5 



/OO /'OO hi 

dW j t/nrjexp -— {/M} 

/ da sin^ a | (r^ + /^) sin^ a + /^ cos^ a \ 



(3.3.7) 



Integrating over a |28| the partition function can be rewritten as. 



Zo 



234^5 ^5^5 

J" 



"" / -n PJT^rt^xp 

-00 Jo {ri + pyi^ 

1 1 






If we scale r^ and / by a factor /5 ^gi, the integral reduces to 



23415 ^^TT^ 

/32 



(irirj 



dll^ 
00 7o (r2 + /2)7/2 



exp 



— r. 



The integral over ri can be done to give 



2^n5 ^-TT 



^^5 



(3.3.8) 



(3.3.9) 



/3 



I" rf//«exp (/Vs) la-^T C-, (/VS)] - 3/^01 (-^, (1^8) 



+3/Vr(^,(/78))-/Vr(^,(/78) 



(3.3.10) 



where r{a,x) is the incomplete Gamma function. 



In large / regime using the asymptotic expression for the incomplete Gamma function ||2^ 
we can write this expression as 

23415 ^f ^5 ^oo^^^^_(4) 



P 



rM-1 



dlfe 



y^ (^_iymg{f +l)/(-2m-3) 



Lm=0 



'r(i +m) _ Jil + rn) ^ ^r(| + m) r(| + m) 



r(i) 



r(|) 



r(|; 



r(i) 



o(|(/V8)r 



-M 



In large / approximation this boils down to 
23^15 g^n^ r^ 



(32 



/oo , 

dll^ (24576/-^^ - 2752512/-^^ + .... 
-oo 



(3.3.11) 



and in the small / regime the above partition function becomes 28 



2^*15 g-ir'' f°° 



dll^Vsr^ exp 



V^ll- 24/^ + 256/^ 



5120 



12 



(-l)"V8/(2+2n) 



n\ 



+ 12E 

I n=0 

and / — i> limit gives 



^ 2"15gT.' .» ^^^, 



(4n2 - 8n - 5)(4n2 - 8n + 3)^ 



(3.3.12) 



P 



^'--'f'<fs-'fy- 



(3.3.13) 



drirj" 



The integrand /* /o°° (^2^2)7/2 exp 



(?)^? 



has /~^ dependence for large / and f dependence 
for small /. Hence, it converges for both large / and small / and the integral is non-singular 
and independent of f3. So, Zq has a temperature dependence of T~ . 



3.4. Non-leading Interaction Term. 



In the previous section we have seen that the leading order fermionic contribution comes from 
the free fermionic terms. Here we will try to estimate the f3 dependence of the non-leading 
fermionic contributions. 
In terms of the parametrisation in eqns. (|3.3.3D , 



Zf = —^ n ( 4^^ + ^^^1^2^ sin a + -^ (ri^ + ra^ + (xif + {X2Y + f 



9 



VT" 



16 



(3.4.1) 



which, in /3 — > can be written as 



^/ = 9-^' n (^'0(/?V) + 0{p-^g-^) + rO{^-*g-^) + 

= (7-2^ (0(1) + 0{pig-^) + 0(/5l^3) + rO{(3-^gi) + ) (3.4.2) 

where the first term is the leading order fermionic partition function we have discussed in 
subsection (3.2). 

So we can see that in the partition function at high temperature the contribution of the zero 
modes (bosonic) is dominant. We have earher argued that the higher modes of the bosonic 
fields will also contribute in non-leading terms. 

3.5. Mean-square Separation of the D-0 branes. 

As we are working in Euclidean metric and since for zero mode calculation we have 
Lorentz symmetry, we can identify / as one of the spatial components and hence as the 
separation between two D-0 branes. 
Now, we try to see the temperature dependence of the expectation value of Z^" 

in = ^^^ (3.5.1) 

i.e. 

foo ^//8/(2n) TOO dr^rj \ //2 \ o" 

roc dll^ C j0^ exp [- (f ) rf] 
As / is the separation of two D-0 branes, we get the mean square separation of two D-0 
branes from this, by putting n = 1 and doing the integral numerically, and restoring Ig, we 

{(a)--(i)' 

If we assume high temperature expression has a finite radius of convergence, we can conclude 
that the mean square separation is finite for finite temperature. This implies that there is a 
confining potential that binds the D-0 branes. As argued earlier the scaling argument that 



gives the P and g dependence in |3.5.3| is valid for all N . So we can conclude that (/^) ~ ^/¥ 



for all N. 
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3.6. Effective Potential. 



For high temperature we have evaluated the partition function both for large and small 
/ (eqn. |3.3.ri| , p.3.130 . Up to leading order the effective potential between two D-0 branes is 



proportional to — log / and log / for small and large /. We can see that the potential increases 
at both I ends, though we can not clearly see the nature of the potential in the intermediate 
region but we can conclude that the potential is a confining potential and binds the D-0 
branes. 

4. CONCLUSION 

In this paper we have attempted to solve the N = 2 matrix model in the high temperature 
limit . The leading nontrivial term of the partition function has been calculated exactly (eqn. 
|3.3.10|) . The non- leading terms can also be systematically calculated although we haven't 
attempted to work them out in this paper. From a scaling argument we have also determined 
the f3 and g dependence of the leading term for any A^. This complements the work of [Q , 
where the one loop partition function was calculated with the entire /3 dependence. We 
find that (/^) oc ^/| (eqn. ^.5.3| ) (true for any A^), the finiteness of which shows that there 
must be a potential between D-0 branes that binds them. In 0J^ also a logarithmic and 
attractive potential were found. The present calculation being exact in Qs is valid at all 
distances. Thus unlike in [§,§, the (finite temperature) logarithmic potential found here is 
attractive at long distances and repulsive at short distances so it has a minimum at non- zero 
separation. Similar issue for low temperature has been discussed in [^]. In [^ it was found 



that at high temperatures, the configuration with all the D-0-branes clustered at the origin 
i.e. with the zero separation, had lower free energy than the one where they were spread 
out. However, that was a large A^ calculation and also restricted to one loop. It is therefore 
possible that more exact calculation will resolve this issue. 

As mentioned in the introduction, describing completely the dynamics of two D-0 branes 
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in M-theory would require the infinite N model. Whether some high temperature expansion 
of that model within the -^ approximation scheme can be attempted is an open question. 

APPENDIX A: COMPARITIVE f3 DEPENDENCE OF THE TERMS IN ACTION 

We are interested in investigating the comparitive /3 dependence of the terms in the 
action in eqn. ( p.0.6|) . For convenience, in this part we will suppress the isospin and vector 
indices. The action takes the form 

1 °° a„ /„ °° 

S* = - < — 2^ — X„X_„ + 2^ —XiXmXn — C 2^ XnXmXiXp > (A.l) 

^ I n— — oo; /^ ^+m+n^O '"^ n,m,Lp— — oo, I 

I, n^O n+m+l+p=0 ) 

where an,bn,c d and /„ are constants and are given by a„ = vr^n^ , c = |, and fn = ^ 
Now, when we expand the sum over n, m, I and p in last two terms, we will get terms with 
all of these indices being 0, with two of the indices being and with one of them being 0, 
so the action can be written in the form (taking one of each type, as the terms of the same 
type have same /? dependence). 

S = --J- 2J (^nX^ -^0 zJ -^O^n + ~ zJ fnXoX_nXn 

r^y n— — 00; J J n^O ^ n^O 

/, XoXmXiXp y^ XnXmXiXp 

Q — Q — 

m+l+p=0 n+m+l+p=0 

Now, we set 

{Xn)nj^o = \J gl3X'^ 
Under this transformation, the action reduces to 

S = (-^Xo4 + \-a^Y: K' - cP' E K'XI + E /„XoX_„X„ 

y 9 ) y „^o n^O n^O 

-c^'i^ E X'^XlX[X,-cP'g J2 KXLX'iX'l (A.2) 

n+m+l=0 n+m+l+p=0 J 
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Again, we can neglect /JaX^X^X/Xo and P^X'^X'^X'^Xp terms also, as they have higher 
order (3 dependence ( | and 3 respectively). To justify this step, note that since the potential 
is always greater than zero, it is equivalent to an integral of the form 

Z = dx exp -ax^ + --^x - (3'^cx^ - cg(3^x'^ = / rfxe""^^ exp ^--^x - (3'^cx^ - cg(3^x^ 

Being a convergent integral, as /3 — *> 0, we can use Taylor expansion and write as 

Z = I e-"-" I 1 - ^x + (3\x' - ^x + ....] dx (A.3) 



V9 4y^ 

Neglecting the higher order f3 terms, we can write the partition function for action ( [A. 21 ), 



Z = JdXoexpl-^X^\ n JdX^exp U-a^ + /?/„Xo - c/3^X^)^xA 

In order to justify dropping the X^X^ and XqX^ terms, we can proceed as follows. After 
doing the X„ integral the partition function reduces to 



exp 



{-f^ol 



As here 4(a„ — l)c — /^ is positive, (a„ — 1) — fn^o + cXq is also positive , i.e. 
[a„-/„Xo + cX2]>l. So 



{an-fnPXo + cP^X^)-2 [ 9 

and the integral is finite. So in /5 ^ limit the leading order term can be given by 

/(iXofln ^ exp I — — Xq > and the corrections will vanish at the positive power of P . So, 

in j3 —>■ limit we can ignore those, which is equivalent to ignoring the fnPXoXnXm and 

c/?^XqX^ terms in action so it can be written as 

oo 

S = - E Xl-cpX', (A.6) 

'^ n=—oo;nj^Q 



which is also eqn. (|3.3.1| ). 
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